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Abstract 



Out of the fifteen 3x3 neutrino mass matrices with two texture zeros seven are compatible 
with the neutrino oscillation data. While 2 of them correspond to hierarchical neutrino masses 
and 1 to an inverted hierarchy, the remaining 4 correspond to degenerate masses. Moreover only 
the first 3 of the 7 mass matrices are compatible with the maximal mixing angle of atmospheric 
neutrino and hence favoured by data. We give compact expressions for mass matrices in terms 
of mass eigenvalues and study phenomenological implications for the 7 cases. Similarity of the 
textures of the neutrino, charged-lepton mass matrices with those of quark mass matrices is also 
discussed. 
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I. Introduction 

There has been a long standing interest in the texture zeros of the 3x3 quark mass matrices 
as a possible source of the observed hierarchies in their masses and mixing angles. In particular 
the up and down quark mass matrices with 2 or 3 texture zeros arc known to successfully relate 
the ratio of quark masses to the mixing angles [1] . In the context of neutrino physics also there is 
a long history of invoking texture zeros for describing the 4 x 4 [2] as well as the standard 3x3 
neutrino [3,4,5] mass matrices. Recently Frampton, Glashow and Marfatia [6] have systematically 
compared the predictions of all the symmetric 3x3 neutrino mass matrices with two or more 
independent texture zeros with the neutrino mass and mixing parameters as derived from the 
oscillation data. A symmetric 3x3 mass matrix has in general 6 independent elements. They find 
no neutrino mass matrix with three or more texture zeros, which is compatible with the neutrino 
oscillation data. Moreover they find that 7 out of the 15 independent neutrino mass matrices with 
two textutrc zeros (*'C2 = 15) are compatible with the oscillation data. The present work carries 
this investigation one step further, giving explicit expressions for the mass matrices in terms of 
the three mass eigenvalues. We also give numerical estimates of these mass eigen values for the 7 
cases, which arc compatible with the neutrino oscillation data. 

As usual the neutrino mass matrices shall be written in their flavour basis, which corresponds 
to the mass basis of the charged leptons. For simplicity we shall consider real mass matrices. 
One can easily check that the main results arc unaffected by the introduction of CP violating 
phases. On the other hand the constraints on these phases arising from the texture zeros have 
been recently discussed in [7]. A 3 x 3 real symmetric mass matrix with two texture zeros has four 
independent parameters. Three of them can be determined in terms of the three mass eigenvalues 
using the invariance of trace and determinant, i.e. 

Tr. = TOi + m2 + ws, Det. M. = TOim2m3, 

Det. Ai X Tr. Jvl^^ = mim2 + m2'mz + ra^rai. (1) 

But we need one experimental input to determine the remaining parameter. We shall assume the 
maximal mixing angle for the atmospheric neutrino oscillation [8] to provide this input, i.e. 

ii = tan 6*23 = 1. (2) 

This represents by far the most robust result of the neutrino oscillation experiments so far. The 
most favoured value of this mixing parameter has remained 1 over the years, while the error bar 
has shrunk steadily. Therefore we treat this as a manifestation of an underlying symmetry of 
the neutrino mass matrix. It is conceivable of course that the physical value of this parameter 
may not be exactly 1. Nonetheless it is fair to assume that this point is smoothly connected to 
the symmetry limit. Therefore we expect the phenomenologically favoured mass matrix to be 
compatible with the maximal mixing angle of eq. (2). We shall see in section III below that 4 out 
of the 7 experimentally allowed mass matrices with two texture zeros are incompatible with this 
requirement. 

The relation (1) was successfully used in determining the structure of the up- and down- 
quark mass matrices in terms of the mass eigenvalues for two and three texture zeroes [9] . The 
supplemental relation in this case came from the assumption that the three angles were very small 
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(mass hierarchy) so that the triangular matrix technique could be implemented [9], [10]. We 
have briefly discussed this in section IV and compared their textures with those of the neutrino, 
charged-lepton system. 

II. Experimental Constraints 
We shall impose the following experimental constraints on the neutrino masses and mixing 
angles. 

i) The atmospheric neutrino data implies 

A„ = Im^ -m2 2| = (1.7- 4)10-3 eV^, (3) 

tan^ 6123 = 4 = 0.9-1, (4) 

at 90% CL [8]. 

ii) The solar neutrino data admits LMS (LOW) solution at 90 (99.7%) CL with [11] 

A, = \ml - mil ~ 6 X 10"^ (1 x lO"'^) eV^ (5) 

^3 = tan=^ 6*12 ~ 0.4 (0.6). (6) 

iii) The CHOOZ and Paoloverde atomic reactor experiments give the 90% CL limit [12] 

S2 = sin ^13 < 0.16. (7) 

Thus we shall require the ratio of the solar and atmospheric neutrino mass differences to satisfy 
the inequality 

A,/Aa = \ml-ml\/\ml-ml2\ < 2 x lO'^. (8) 

III. Neutrino Mass Matrices 

The mass matrix can be expressed in terms of the above masses and mixing angles as 

/mi \ 

M = U I m.2 f/"^, (9) 

V ma/ 

where 

C2C3 C2S3 S2 \ 

U = I -C1S3 - S1S2C3 C1C3 - S1S2S3 S1C2 I . (10) 
S1S3 - C1S2C3 -S1C3-C1S2S3 C1C2/ 
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Using the maximal mixing constraint (2) wc get the following expressions for the matrix elements 
of M., which are valid upto the first order terms in S2- 

A4ii = mic\ + m2s\, 

Mi2 = -{mi - m2)s3C3/\/2 + m3S2/V2 - {micl + m2sl)s2/V2, 
Mi3 = {mi-m2)s3C3/V2 + m3S2/V2- {micl + m2sl)s2/V2, 

M22 = (miSg +TO2C3)/2 + TO3/2 + (mi - m2)S2S3C3, 

M33 = (miSg + m2C3)/2 + 7713/2 - (mi - m2)s2S3C3, 

M23 = -(miS3 +m2C3)/2 + m.3/2. (11) 

From these expressions one can study the implications of setting different pairs of matrix elements 
to zero, as we see below. 

(A) Hierarchical Solutions 

1) A^ii = 0, ;Wi2 = 0: They imply 



= -mi/m2, S2 = V-mim2/m3. (12) 

Combining these with eqs. (6) and (7) we get |mi| < |m2| <C m3, i.e. hierarchical masses. 
Prom eqs. (3) and (5) we get Ims] ~ .05, |m2| =i .009, |mi| ~ .004 eV for the LMA solution; 
while the mi and m2 values are each suppressed by a little over an order of magnitude for 
the LOW solution. In the former case we get S2 — 0.12, i.e. close to the CHOOZ limit (7); 
while it is over an order of magnitude smaller in the latter case. Note that for this solution 
the Ve Majorana mass, 

Mn = Y.U^,m,, (13) 

is zero. Thus it predicts no 0i//3/3 signal even at the level of .001 eV, which corresponds to 
the highest level of sensitivity expected at the proposed GENIUS experiment [13] . 

Finally substituting (12) in (11) gives the explicit form of the mass matrix to first order in 
S2, i-e. 

/ y/-2mim2\ 

_ I Q mi+m2+m3 m3-mi-m2 j {^^) 

VV-2miTO2 "^-^-""^ nii+A+imJ 

One can also derive this directly from eqs. (1) and (2). A special case of this mass matrix 
corresponding to bimaximal mixing, ts = —1711/171,2 = 1, has been obtained in a U(l) gauge 
extension of the standard model, corresponding to the gauge charge — [14] . However we 
could find no simple dynamical model for the general case of mi ^ —m2- Some explorations 
in this direction can be found in refs. [5] and [15]. 

2) Mil — 0, Mi3 — 0: The phenomenological predictions for this case are identical to the 
previous one except for the change of sign of S2- The mass matrix is simply obtained from 
(14) by interchanging the 12 and 13 elements. 
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(B) Degenerate Solutions 



1) M.22 = 0) -^13 = 0: Substituting these in eq. (11) imply 

m\ = m\ + 4:S2ml/t3, m\ = m\ — 4s2TO3t3, (15) 

i.e. ra\, m§, are nearly degenerate and occur in that order. Thus \m^ — ra\\ > [m^ — m2|, 
in gross contradiction with the observed inequality of the solar and atmospheric neutrino 
masses (eq. 8). 

As shown in [6] , this mass matrix is compatible with the experimental constraint of eq. (8) 
away from the maximal mixing points (2), i.e. for ti ^ 1. In this case one gets 

ml = ml[4 + 2s24{ti + l/h)/t3], 

ml = ml[4-2s2tl{ti + l/ti)t3]. (16) 
Compatibiity with the experimental constraint of eq. (8) can be achieved for 

S2 ~ 0.5 X 10-2(1 - t1)/{t3 + l/ts) < 4 X 10-^ (17) 

corresponding to the > 0.9 range of eq. (4). Thus the S2 angle in this case is too small to 
be observed at the future long base line experiments. On the other hand the Majorana 
mass 

Mn ^ 4^1^^^ > 0.10 eV. (18) 

This is not too far below the present experimental upper limit of 0.2 eV [16], and will surely be 
measurable at GENIUS [13]. Note that this represents the common mass scale of degenerate 
neutrinos and it is cosmologically significant. Nonetheless we consider the incompatibility of 
this solution with the maximal mixing region {ti = 1), favoured by the atmospheric neutrino 
data, to be a serious drawback of this model. The expressions for the matrix elements are 
rather long for ti ^ 1. Therefore we are not displaying the explicit form of the mass matrix. 

2) M.S3 = 0, M.12 = 0: Substituting these in eq. (11) imply 

ml = TOg — 4:S2ml/t3, m\ = ml + 4s2mlt3, (19) 

i.e. nearly degenerate , m| , m| with |m| — wfl > |to| — as in the previous case. In fact 
the magnitudes of S2 and t3 are the same as above. The results for ti ^ 1 are similar to the 
previous case with ti replaced by — 1/ti. Consequently the predicted S2 and Vg Majorana 

mass are very close to those of eqs. (17) and (18). 

3) A422 — 0, Aii2 = 0: The results are practically the same as in 2). 

4) A433 = 0, A4i3 = 0: The results are practically the same as in 1). 

Matrices Bi through B4 in terms of mass eigenvalues are listed in Appendix I. 
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(C) Inverted Hierarchy 

M.22 = 0, A^33 = 0: Substituting these in eq. (11) gives 

S2 = 0, tl = -{1712 + ma) /{mi + ma). (20) 

The first cquaUty impHcs of no observable CP violation in the neutrino sector. The consistency of 
the second with eqs. (6) and (8) implies an inverted hierarchy of masses, i.e. 

-m2 mi ~ (2 - 4)m3. (21) 

Substituting this in eq. (3) gives — m2 ~ mi ~ .05 eV and ms ~ .02 eV. The predicted Ve 
Majorana mass is 

1 _ +2 

Mil — ^\/a;~ (.025-.013) eV. (22) 

This is an order of magnitude below the present upper limit of 0.2 eV [15], but will be measurable 
at GENIUS [13]. Substituting (20) in (11) gives the explicit form of the mass matrix 



mi + m2 + TOa -y/-{mi + ma){m2 + m3)/2 y^-(mi + mz){m2 + m3)/2^ 

M = I -a /- (mi + m3,){m2 + m3)/2 ms 

-\/-(mi + ma){m2 + m3)/2 ms 

(23) 

(D) Disallowed Cases 

Finally let us briefly indicate why the remaining 8 mass matrices are experimentally disallowed. 

1) Mi2 = 0, A^is = 0: A^i2 - Mia = implies 

(mi - m2) sin 26*12 = 0, (24) 

i.e. no solar neutrino oscillation. 

2) A^i2 or Mi3 = 0, M23 = 0: V27Wi2(i3) - 2s27W23 = implies 

I tan 261121 = = 2\s2\ < 0.32, (25) 

in gross disagreement with the solar neutrino result (6). 

3) Mil = 0, M22 or M33 or M23 — 0: In each case one gets (neglecting 0{s2) terms) 

jm^ - m^l = (1 - 4)ml, \ml - m^j = {24 - 4)ml, (26) 
in conflict with eqs. (6) and (8). 

4) M22 or M33 = 0, M23 = 0: 22(33) + M23 = implies 

ms = - (mi - m2)s2S3C3, i-e. mi ~ -m2 » ms (27) 
using eqs. (7) and (8). Substituting this in M23 = implies ~ 1, in conflict with eq. (6). 
Matrices Di through Dg in terms of mass eigenvalues are listed in Appendix I. 
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IV. Comparison with Quark Mass Matrices 



As stated earlier, equation (1) was used in obtaining the structure of the up- and down- 
quark mass matrices, U and D, with two and three texture zeroes, in terms of the mass eigenvalues 
[9]. In determining their matrix elements it was further assumed that the mixing angles were 
very small so as to be compatible with mass hierarchy [9], [10]. Among the available choices, the 
structures that fit the experiment the best, with two texture zeroes, were [9] 



U = 









with mi 
and 



m2 



ma = mt 



s/-mim2 



(28) 







D = I v'— mim2 




V'— mim2 

m2 ^mims 
s/mimz ms 



(29) 



with mi = TUd, m2 = m^, ms = m^. 

These matrices together give the correct CKM angles [9] 



It is interesting that the texture of the above U - matrix is similar to the best choice for the 
neutrino mass matrix {M.) for S2 > , and hierarchical mass given by (14). One may conjecture 
that in that case the texture of the charged-lepton mass matrix (L) should be the same as the 
£>-matrix. Indeed, since the hierarchy of masses in L is similar to D, the entire structure of the 
two matrices may be similar. 

In other words, the best choice for the neutrino and charged-lepton mass matrices would 
be (14) and 





L = I ^-mem^ 




rrin 







(30) 



We realize that the neutrino mass matrices were written above in the basis in which L is 
diagonal whereas the above L is not. One could diagonalize L by rotating it through appropriate 
angles. The new A4 matrix will then be different from matrix (14) by terms proportional to the 
rotation angles which, however, are very small compared to their neutrino counterparts since they 
arc proportional to the ratios of the charged lepton masses which arc, indeed, very small. Even 
for which is proportional to S2 , and, therefore, small, the new AI13 will not be drastically 

different since the relevent rotation angle here is proportional to L13 which is actually zero. Thus 
we expect the new A4 not to be very different from (14), inspite of the fact that (30) is not 
diagonal. 
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V. Conclusion 



Only a limited number of texture patterns for the neutrino mass matrices are consistent with 

experiment, particularly if one takes the experimentally favoured value of ti = 1. These are given 
by ^1,^2, and C . In particular, for hierarchical masses and S2 > there is only one possible 
matrix Ai given by (14). We also note that this particular matrix has a texture which is similar 
to the up-quark matrix, that is found compatible with experiments, which leads to the possibility 
that the quark and lepton mass matrices may have similar texture patterns. 

We thank Anjan Joshipura and Ernest Ma for discussions. This work was supported in part 
by the U.S. Department of Energy under Grant No. DE-FG03-94ER40837. 
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Appendix I 



For the sake of completeness we compile below the neutrino mass matrices for ti = I for the 
cases B and D, although, as discussed earlier, they are in conflict with the neutrino oscillation data. 



( mi + TO2 + ma -\/-2(mi + m,3)(m2 + ma) \ 

^-2(mi + m3)(m2 + ma) ma j 

ma / 
(Note: extra texture zero develops since 0(s2) are neglected.) 



(mi + m2 + ma ^y-2{ml + ma)(m2 + ma) 
ma 
^-2(mi + ma)(m2 + ma) ma 
(Note: extra texture zero develops since 0(s|) are neglected.) 



^i) 



D2) 



' mi 

n m3+"i2 
" 2 

ri m3— 7rt2 
" 2 



yn3 — m2 
2 

m3+m2 
2 



mi 



m2 




ma 



-\/2-\/ma(mi +m2 - 2ma) - (mi 



m2 



m§) 



mim2 — y».3 
mi+m2— 2m3 





—\f2^ m-i{m\ +1712 - 2m3) - (mim2 - mD^ 


2m3(mi+m2— 2m3) — (^17712— 
mi+m2 — 2m3 



mi + 777.2 ~ 

Ds) I A/2-\/m3(mi + m2 - 27773) (mim2 - m|) 





V2^/ms{ml + m2 - 27713) - (mim2 - ml) 

2m3 (mi +m2 — 2m3 ) — (mi m2 — ^ris ) 
mi+m2— 2m3 








m\m2—Tn\ 
mi+m2 — 2m; 



Da 



( 

— 2m2mi +m3 (mi +m2+m3) 

2v^— 2mi7n2 
2m2mi+m3(mi+m2+m3) 

\ 2^/— 2mim2 



-2m2mi+m3(mi+m2+m3) 
2\/— 2mim2 



i(-mi - m2 + ma) 



2m2TOl+TO3(TOl+TO2+TO3) \ 

2v'— 2mim2 

i(-m,i 



mi 



— m,2 + m,3j 
m2 + ma 



( 

-2m2mi — m3(mi+m2+m3) 

2V — 2mim2 
2m2mi—m3{mi+m2+m3) 

\ 2v'-2mi 



-2m2mi — m3(mi+m2+m3) 2to2TOi — m3(mi +m2+m3) \ 
2y/ —2mim2 2-^— 2mim2 



mi + m2 + ma 
i(-mi - m2 + ma) 



5 (-mi - m2 + ma) 








/ -mim2 _ /- 



mi + m2 

-TOim2 Q 



-mim2 





mi + m2 
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^ I ™ _ rn TOim2-TO3 2to3 (toi +TO2 -2m3) - (mim2 -mg ) \ 

ffbi ~\~ f f t'2 ' ' *'3 / / „ — , I — 

V^y m3(mi+m2 — 27713) — (77117712—7713) V^y 7713(7711+7712— 27713) — (77117712 — 







V^-^ m3(mi+m2 — 27713) — (mim2—™l) 

2m3(mi+m2 —2m;i} — (inim2 —"13) q 
\ \/2-\/ m3(Tni+m2 — 2m3) — (mim2— m|) y 



^ ^ -j- 777, 777, — 2m3(mi+m2— 2m3) + (mim2— 7713) m^— ?ni7n2 ^ 

V2-\/ "i3('"l+™2 — 2m3) — (mim2— m^) >/2-\/"i3 ("^l +"^2 — 2m3) — (mim2 — rrig) 





- 2m3(mi+r?T,2 —27713 )4-(y>T'i?Ti2—y^3) 



V2-yrn3(r7ii+r?i2 — 2r?i3) — (r7iir?i2— r^l) 

777,3- TTJ-i 7772 



2m3 



y \/2y' m3(mi+m2 — 2m3) — (mim2— mg) y 
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